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Abstract 



A remarkable orbital quadrupole magnetic resonance, so-called twist mode, is 
predicted in alkali metal clusters where it is represented by I T = 2~ low-energy 
excitations of valence electrons with strong M2 transitions to the ground state. 
We treat the twist by both macroscopic and microscopic ways. In the latter 
case, the shell structure of clusters is fully exploited, which is crucial for the 
considered size region (8 < N e < 1314). The energy- weighted sum rule is 
derived for the pseudo-Hamiltonian. In medium and heavy spherical clusters 
the twist dominates over its spin-dipole counterpart and becomes the most 
strong multipole magnetic mode. 
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Orbital magnetism in atomic clusters is a subject of special interest. Clusters may contain 
many atoms and, therefore, single-particle orbital moments of valence electrons can reach 
very big values. This results in huge orbital effects, for example, in strong orbital magnetic 
resonances. These resonances are of a general character and exist in different finite Fermi 
systems (nuclei, atomic clusters, etc). Ml scissor mode in deformed systems pi]- H] and M2 



twist mode in systems of arbitrary shape ||- |L0| are most famous examples. They have 
been observed in atomic nuclei (see, e.g., 0,[J) but not yet in clusters. The scissors mode 
has been already predicted in clusters ||,[4| . In the present paper we will discuss properties 
of the twist mode in this system. 

By definition, the twist is the quadrupole torsional vibration mode of an elastic globe || . 
It is generated by the operator T = e~ tazlz = e au ' v with the velocity field u = (yz, —xz, 0) 
|5|||]. The mode is viewed macroscopically (see Fig. 1) as small-amplitude rotation-like oscil- 
lations of different layers of a system against each other with a rotational angle proportional 
to z (projection to the axis of rotation). The restoring force of the twist is determined by the 
quadrupole distortions of the Fermi surface in the momentum space. So, the twist represents 
transverse magnetic quadrupole oscillations of an elastic medium, provided by variations of 
the kinetic-energy density. The twist is a general feature of any 3-dimensional finite Fermi 
system which demonstrates an elastic behavior. Atomic nuclei ||- @ and clusters [10 
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most typical examples. Unlike the Ml scissor mode which has the similar quantum origin 
but can exist only in deformed systems, the twist manifests itself in Fermi systems of any 
shape, spherical and deformed. 

Expressions for the twist energy and M2 strength, obtained within the elastodynamical 
models [0,R| for atomic nuclei, can be reformulated for clusters as 



17eV A 2 rJ 2 N- 1/3 , B{M2) = 0.52r 2 s N^ 2 b (0.1) 



where B(M2) is probability of M2 transition from the ground state to 2~ twist state, r s 
is the Wigner-Seitz radius, and N e is the number of valence electrons. We will show that 
expressions ( |0.1| ) provide good qualitative agreement with microscopic results. However, 
elastodynamical models do not access the shell structure of nanoparticles and so cannot 
clarify the microscopic origin of the twist. Moreover, these models are questionable for 
small systems. In this paper we will present, for the first time, a microscopic analysis of the 
twist, fully embracing shell effects. Both small and heavy clusters will be covered. 

The operator for M2 transition F(M2,/i) = fi b y/T0r[g s {Y 1 s} 2 ^ + ^gi{Y 1 l} 2 ^}, is a 
sum of spin and orbital components with g s = 2 and gi — 1. The external field generating 
the twist, zl z oc r(Y w l z ), is a part of the orbital component with fi = 0. So, it is natural to 
consider the twist as a part of the orbital M2 resonance. Since both are of a similar nature, 
we will call the whole orbital M2 resonance also as a twist. 

In spherical systems, only the spin-dipole channel delivers the residual interactions for 
2~ excitations. Investigations 0,0 have shown that twist in spherical atomic nuclei only 



weakly depends on this residual interaction and the remaining dependence is mainly caused 
by the spin-orbit coupling. The influence of the residual interaction should then be even 
smaller for metal clusters where the spin-orbit coupling is negligible. Thus we can deal with a 
simple particle-hole (lp-lh) picture. Also, we will confine our consideration to the spherical 
jellium model which is compulsory for analysis of large systems. These two approximations 
are appropriate for the present purpose of a first survey of the twist mode. At the side of 
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the mean field Hamiltonian, we take into account local as well as non-local effects which are 
caused by the ionic pseudo-potential. For Na and K clusters we calculate the Kohn-Sham 
single-particle scheme within the approach [13j which properly treats the local effects. For 
Li clusters, the pseudo-Hamiltonian [|T5]| 



Ti 2 - - -* 

H = -—- V ■ (1 + a(r))v + L ■ (3{r)L + u{r) + W(r) (0.2) 

— 1 1 tf> 



with the parametrization [TBI is used. The local ionic contribution is carried in u{f) and 



non-local effects lead to the effective mass m*(r) = m e /(l + a(r)) and the orbital term 
oc (3{r). Coulomb and exchange-correlation potentials are represented by W(r). 

Results of our studies are exhibited in Figs. 2-4 and Table 1. In Fig. 2 the orbital M2 
strength, B(M2) = \ < i\F(M2)\0 > | 2 , is presented for light, medium and heavy spherical 
Na and K clusters (i represents (lp—lh) 2 - states). Only the orbital part of the M2 transition 
operator is used. Fig. 2 shows that dominant twist strength is concentrated in the one lp-lh 
peak with the lowest energy. This relation persists independent of N because both the twist 
mode and other relevant lp-lh excitations follow a trend oc iV" 1 / 3 (see also Fig. 4). The 
degree of concentration, however, changes with N. The twist peak exhausts 100%, 80%, and 
60% of the total M2 strength in Nag , Ndg 3 , and Naf 315 , respectively. It corresponds to the 
nodeless branch of the twist mode (see Fig. 1, left side). The nature of the peak is clarified 
in Fig. 3. It represents n, I — > n, I + 1 transition between single-particle levels with the node 
number n — 1 and maximal orbital moments I. The levels belong to the last occupied and 
first empty shells. As a result, the twist can serve as a valuable source of information about 
i) single-particle levels with maximal orbital momenta near the Fermi surface and ii) the 
energy gap, AE s h, between the Fermi and next empty shells. 

In large clusters weaker peaks also contribute to the twist resonance. As a rule, they 
represent n, I — > n, I + 1 transitions with n = 2, 3, ... and lower orbital moments. Following 
our analysis of the velocity field, these peaks also mainly contribute to the nodeless twist 
branch (Fig. 1, left side). Their contribution grows with increasing cluster size. They come 
energetically closer to the dominant peak (see Fig. 2) such that they may not be easily 
distinguishable experimentally. Other twist branches (see, e.g., the right side of Fig. 1 ) 
carry only a small fraction of the total strength and lie at higher energies. 

Fig. 2 demonstrates similarity of the twist in Na and K clusters. Much similar pattern 
appears for Li. However, there are also distinctive differences in twist energies and strength. 
These values are compared in Fig. 4. K, Na and Li are distinguished by Wigner-Seitz radii 
(in atomic units, r s (K) = 5, r s (Na) = 3.96 and r s (Li) = 3.25, respectively) and by different 
influence of the ionic structure. Fig. 4 shows the following trends: i) The denser the metal, 
the higher the twist energy. This can be explained by the fact that the smaller the r s , the 
deeper the corresponding single-particle potential |I4] |. In Li clusters the potential is most 



deep and, therefore, has most large energy gap AE sh between neighboring quantum shells. 
Being close to AE s h, the twist energy should increase from K to Li. ii) In all the size region 
the twist mode resides far below the Mie dipole plasmon and, at the same time, it stays 
still far above typical ionic vibrations (u ~ 50 meV), which both helps for an experimental 
discrimination. The heavier the cluster, the smaller the energy gap AE s h and, therefore, 
the twist energy. Fig. 4 demonstrates the energy fall oc iV" 1 / 3 . It is worth noting that such 
dependence of twist energy takes place in both atomic clusters and nuclei, iii) The denser 
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the metal, the smaller the resonance strength. As is shown below, the energy-weighted sum 
rule Si keeps the same value for K, Na and Li clusters of a given size (if one neglects the 
ionic structure effects). So, the increase in the excitation energy has to result in the decrease 
of the strength. In Li clusters the strength is additionally suppressed due to the effective 
mass, m*/m e ~ 1.2. iv) B(M2) grows as N% (see upper part of Fig. 4). 

All the trends discussed above are supported by the elastodynamical results ( p.l|) . How- 
ever, the quantitative agreement is less perfect. Eqs. ( p.l| ) considerably overestimate both 
twist energy and strength (up to 80% and 30% in light and heavy clusters, respectively). 
That is not so surprising because we find that the twist mode is dominated by the shell 
structure while the collective model averages over shells. 

The attractive feature of the orbital M2 mode is that its total strength can be estimated 
in a simple fashion by thwe energy-weighted sum rule as 

S = J2°°i\ < i\F(M2)\Q > | 2 (0.3) 

i 

l 2 



where 



- J2 (0\[F(M2 fl ),[H ,F(M2 f i)]]\0) = S 3 + S l 

A fi=-2 



7^ 2 ""^ r 
S s = « E( 2 ^ + 1) / (1 + « + 8r 2 (3)p nth dr, (0.4) 

Si = -P- E(2/, 3 + 3/ 2 + k) / (1 + a + \r*P)p nili dr 
6Trm e , J 5 

are contribution of the spin and orbital parts of the transition operator F(M2). In Eq. ( p.3|) 
the sum runs over all (lp — lh) 2 - states. In Eqs. ( p.4|) , a(r) and /3(r) are the functions 
responsible for the nonlocal effects in the pseudo-Hamiltonian (p.2|), p ni i t = (r_R n ^. (r)) 2 , 
R n i{r) is radial wave function of the single-particle state nl. In Eqs.( p.4|) , the sum runs over 
all occupied single-particle levels. The spin-orbit coupling is neglected in the calculations, 
but it would anyway not contribute directly to Si. Without the spin-orbit coupling the 
spin-dipole residual two-body interaction also does not contribute to Si. The calculations 
show that even for Li clusters one may safely neglect the spin-orbit correction (3 because 
it contributes less than 0.1%. The influence of the correction a is strong in Li: we have 



m* e /m e ~ 1.2 and so this correction decreases S\ by a factor of ~ 4/5. 

If one neglects the non-local corrections (which is justified for Na and K clusters [P"5|,|T7|] ) , 
then the radial integrals in Eqs. (|0.4| ) are just 1 and both S s and Si become very simple and 
even model independent (the similar sum rules have been derived for the twist in atomic nuclei 
0] and spin-dipole excitations in clusters |I5| ). This is evident for S s where 2X^ (2^ + 1) = 
N e . As for Si, one should mention that different models predict, as a rule, the same sequence 
of occupied levels in spherical clusters, at least for light and medium sizes. Finally, it is worth 
to emphasize that: i) the calculation of Si is extremely simple: it is enough to know the 
orbital moments I of occupied single-particle levels; ii) both S s and Si values are equal for 
clusters of a given size but of different metals (K, Na), which agrees with Eqs. ( p.l| ). 



4 



The twist part, fib \J '80/ '27 'rY w l z , of the transition operator F(M2) gives exactly 4/9 of 
the complete Sj values. Following Fig. 1, the twist represents the oscillations in equatorial 
planes. The total orbital M2 resonance takes also into account the meridian oscillations. 

Results presented in Table 1 show that the orbital M2 energy-weighted strength domi- 
nates over the spin one already in clusters with N e = 40. Starting with N e = 92, the orbital 
contribution becomes overwhelming and, already for N e = 440, demonstrates a huge value 
of 2 • 10 5 fi b l A 2 eV . Since the long-wave Ml response, both spin and orbital, is forbidden 
in spherical clusters (indeed, both spin and orbital lp-lh matrix elements of Ml transi- 
tion are proportional to the radial integral / it!„ 1 ; 1 (r)it! n2 ; 2 r 2 dr = S ni i uri2 i 2 which is zero in 
the nondiagonal case due to orthonormalization condition), the twist starting with medium 
sizes becomes the strongest multipole magnetic mode. This fact emphasizes its fundamental 
character. 

Our calculations indicate that twist mode cannot be detected in photoabsorption spectra 
since it is masked by low-energy El excitations. These excitations are much weaker than the 
dipole plasmon but, nevertheless, strong enough to mask the twist. The inelastic scattering 
of polarized optical photons (resonant Raman scattering) seems to be more appropriate to 
observe the twist, though any conclusions about perspectives of these reactions still requires 
a careful analysis of the competition between El, E2 and M2 modes. The Raman scattering 
can separate electric and magnetic modes due to the polarization selection rules. Clusters 
with about 10 4 atoms seem to be optimal. In such clusters the twist strength reaches 
impressive values and, at the same time, is strongly concentrated at very narrow low-energy 
interval (which remains still well separated from ionic vibrations). Our estimations show that 
in heavy clusters, in spite of a dense spectrum, the twist exists and carries the spectroscopic 
information mentioned above. 

In summary, the M2 orbital resonance and its important part, the twist mode, have 
been investigated in spherical alkali metal clusters. The macroscopic treatment of the twist 
exhibits this mode as a general feature of any finite 3-dimensional Fermi system and provides 
a pertinent description of the basic trends. However, it is not delivering a quantitative 
agreement. The microscopic treatment, including the novel energy- weighted sum rules, 
clarifies the main properties of twist lp-lh M2 response. The resonance is mainly exhausted, 
first of all in clusters of a moderate size, by one lp-lh M2 transition connecting the Fermi and 
next empty quantum shell, namely, their levels with maximal orbital moments. As a result, 
the twist can provide a valuable information about the single-particle scheme. Twist energy 
and strength evolve with a cluster size as lu ~ A" 1 / 3 and B(M2) ~ iV e 2 . In heavy clusters 
an impressive M2 strength can be reached. The twist dominates in the low-energy region 
over its spin-dipole counterpart already in clusters of a moderate size and finally becomes 
the strongest magnetic multipole mode. We hope that the fundamental significance of the 
twist for orbital magnetism in spherical clusters will encourage experimentalists to look for 
proper ways for its observation. 
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FIGURE CAPTIONS. 

Figure 1: Nodeless (left) and one-node (right) branches of twist mode ||. 

Figure 2: The distribution of M2-strength in spherical Na (left) and K (right) clusters 
of different sizes, as indicated. 

Figure 3: Single-particle levels and M2 lp-lh transitions in Nag" and Na^. The Fermi 
levels are marked by the double line. For the main transitions (bold arrows), the contribu- 
tions to the complete strength B(M2) are given. 

Figure 4: The strength normalized by iV| (upper panel) and averaged energy (lower 
panel) of the twist resonance in K, Na and Li clusters. The trends B(M2) oc N% and 
u oc iV" 1 / 3 are distinctive. 
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TABLES 



TABLE I. Orbital energy- weighted sum rule Si (in units ^A 2 eL) and ratio R = Si/S s for Na 
and Li spherical clusters in the energy interval 0-6 eV (95 — 99% of the sum rules are exhausted in 
the interval 0-2 eV). Values R are equal for K, Na and Li clusters of the same size. 



20 
40 
92 
440 



Si 



K, Na 
1.2L10 2 
7.28-10 2 
2.55-10 3 
1.35-10 4 
2.01-10 5 



Li 
1.08-10 2 
6.42-10 2 
2.23-10 3 
1.15-10 4 
1.67-10 5 



R 

K, Na, Li 
0.33 
0.80 
1.40 
3.21 
9.72 
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